Search results for "Spectral method"
showing 10 items of 32 documents
Unsteady Separation and Navier-Stokes Solutions at High Reynolds Numbers
2010
We compute the numerical solutions for Navier-Stokes and Prandtl’s equations in the case of a uniform bidimensional flow past an impulsively started disk. The numerical approx- imation is based on a spectral methods imple- mented in a Grid environment. We investigate the relationship between the phenomena of unsteady separation of the flow and the exponential decay of the Fourier spectrum of the solutions. We show that Prandtl’s solution develops a separation singularity in a finite time. Navier-Stokes solutions are computed over a range of Reynolds numbers from 3000 to 50000. We show that the appearance of large gradients of the pressure in the stream- wise direction, reveals that the visc…
Singularity formation for Prandtl’s equations
2009
Abstract We consider Prandtl’s equations for an impulsively started disk and follow the process of the formation of the singularity in the complex plane using the singularity tracking method. We classify Van Dommelen and Shen’s singularity as a cubic root singularity. We introduce a class of initial data, uniformly bounded in H 1 , which have a dipole singularity in the complex plane. These data lead to a solution blow-up whose time can be made arbitrarily short within the class. This is numerical evidence of the ill-posedness of the Prandtl equations in H 1 . The presence of a small viscosity in the streamwise direction changes the behavior of the singularities. They stabilize at a distanc…
Spectral approach to the scattering map for the semi-classical defocusing Davey–Stewartson II equation
2019
International audience; The inverse scattering approach for the defocusing Davey–Stewartson II equation is given by a system of D-bar equations. We present a numerical approach to semi-classical D-bar problems for real analytic rapidly decreasing potentials. We treat the D-bar problem as a complex linear second order integral equation which is solved with discrete Fourier transforms complemented by a regularization of the singular parts by explicit analytic computation. The resulting algebraic equation is solved either by fixed point iterations or GMRES. Several examples for small values of the semi-classical parameter in the system are discussed.
Discretization estimates for an elliptic control problem
1998
An optimal control problem governed by an elliptic equation written in variational form in an abstract functional framework is considered. The control is subject to restrictions. The optimality conditions are established and the Ritz-Galerkin discretization is introduced. If the error estimate corresponding to the elliptic equation is given as a function like where h is the discretization parameter and is an integer, then the error estimates for the optimal control, for the optimal state and for the optimal value are obtained. These results are applied first for a Two-Point BVP and next for a 2D/3D elliptic problem as state equation. Next a spectral method is used in the discretization proc…
Etude numérique d'équations aux dérivées partielles non linéaires et dispersives
2011
Numerical analysis becomes a powerful resource in the study of partial differential equations (PDEs), allowing to illustrate existing theorems and find conjectures. By using sophisticated methods, questions which seem inaccessible before, like rapid oscillations or blow-up of solutions can be addressed in an approached way. Rapid oscillations in solutions are observed in dispersive PDEs without dissipation where solutions of the corresponding PDEs without dispersion present shocks. To solve numerically these oscillations, the use of efficient methods without using artificial numerical dissipation is necessary, in particular in the study of PDEs in some dimensions, done in this work. As stud…
Fast algorithms for free-space diffraction patterns calculation
1999
Here we present a fast algorithm for Fresnel integral calculation. Some fast algorithms using the fast Fourier transform are analysed and their performance has been checked. These methods are of easy implementation, but are only valid for a specific range of distances. Fast algorithms based on the Fractional Fourier transform allow accurate evaluation of the Fresnel integral from object to Fraunhofer domain in a single step.
Controllability method for acoustic scattering with spectral elements
2007
We formulate the Helmholtz equation as an exact controllability problem for the time-dependent wave equation. The problem is then discretized in time domain with central finite difference scheme and in space domain with spectral elements. This approach leads to high accuracy in spatial discretization. Moreover, the spectral element method results in diagonal mass matrices, which makes the time integration of the wave equation highly efficient. After discretization, the exact controllability problem is reformulated as a least-squares problem, which is solved by the conjugate gradient method. We illustrate the method with some numerical experiments, which demonstrate the significant improveme…
Un modelo no hidrostático global con coordenada vertical basada en altura
2015
Memoria presentada en la Universitat de València para optar grado de de Doctor Esta tesis documenta la investigación que he realizado en modelización atmosférica: se parte de las ecuaciones físicas de la atmósfera y se aplican métodos numéricos eficientes para encontrar una solución a dichas ecuaciones a partir de unas condiciones iniciales dadas. Para este fin, se ha desarrollado un modelo atmosférico cuyas características principales son: espectral en la representación horizontal de los campos, discretización vertical de alto orden de exactitud, y semi-implícito en la integración temporal. Además, el modelo es no hidrostático y tiene una coordenada vertical basada en altura, en vez de la …
Fatigue life under non-Gaussian random loading from various models
2004
Fatigue test results on the 10HNAP steel under constant amplitude and random loading with non-Gaussian probability distribution function, zero mean value and wide-band frequency spectrum have been used to compare the life time estimation of the models proposed by Bannantine, Fatemi–Socie, Socie, Wang–Brown, Morel and Łagoda–Macha. Except the Morel proposal which accumulates damage step by step with a proper ethodology, all the other models use a cycle counting method. The rainflow algorithm is used to extract cycles from random histories of damage parameters in time domain. In the last model, where a strain energy density parameter is employed, additionally spectral method is evaluated for …
Fatigue life calculation by means of the cycle counting and spectral methods under multiaxial random loading
2005
The paper contains a new algorithm for estimation of fatigue life in HCF regime under multiaxial random loading using spectral methods. Loading of Gaussian distribution and narrow- and broad-band frequency spectra were assumed. Various characteristic states of multiaxial loading were considered. The equivalent stress history was determined with use of the failure criteria of multiaxial fatigue based on the critical plane. For determination of the critical plane position, the method of variance was applied. During simulation, the authors compared the results obtained by a spectral method in the frequency domain with those from the rain-flow algorithm in the time domain. The paper also contai…